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The many-anyons wavefunction is constructed via the superposition of all the permutations on 
the direct product of single anyon states and its interchange properties are examined. The phase of 
permutation is not a representation but the word metrie of the permutation group . Amazingly the 
interchange phase yields a finite eapaeity of one quantum state interpolating between Fermion and 
Boson and the mutual exchange phase has no explicit effect on statistics. Finite capacity of quantum 
state is defined as Gentile statisties and it is different from the fractional exclusion statistics. Some 
discussion on the general model is also given. 
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In the past three decades, the topic of anyon [l| and 
fractional statistics is continually of great interest for 
physics coniniunit}j2n3| . Anyon is defined through the 
interchange phase quasi-particles with arbitrary 

interchange phase 0 could exist and lead to an interme¬ 
diate statistics Q interpolating between Boson(^ = 0) 
and Fermion(^ = jr). The widely analyzed model is the 
charge-flux mo del [3, 0 and spin lattice [iqI-[i 2|. although 
the interchange property for many-body wavefunction is 
never carefully examined. Alternatively, in condensed 
matter system the one-particle Hilbert space dimension 
da in general depends on the occupation of state, 
Ada = dapANp, yielding fractional exclusion 

statistics proposed by Haldane- wulHlii. The statisti¬ 
cal parameter QajB might relate to the interchange phase, 
e.g. Qa^ = 6>^^/( 27r) in FQHEjlll and Qa^ = in 

spin chain|lQ[ Il4| with possible verification |1 hi llbj. Yet, 
an universal relation between the interchange phase and 
statistics effect is still lacking. The aim of this paper is 
to construct the many-anyons wavefunction and explore 
its statistics. 

Many-anyons wavefunction-The physical interchange 
of anyons may give any (rational) phase 0, (faCfl = 
e^^(j)\(j)\, specifically in two dimension as a representation 
of the braid group [Hi [11 . Analogy to the Slater deter¬ 
minant form of many-Fermions wavefunction, we could 
build up our many-anyons wavefunction(unnormalized) 
by summing over all the permutations on the direct prod¬ 
uct of single particle states: 

^ ( 1 ) 

P 

here <Fo is the direct product of the orthonormal sin¬ 
gle anyon states (fa with a preset order, let’s say 4>o = 
(fa® fb ’ ’ ’ (omit for later convenience). P belongs 
to the permutation group Sn- P^o is a permutation of 
the ordered anyons and its coefficient is a phase 
to guarantee their equal probability. Once we select the 
phase of 4>o as zero, then we will expect the phase of P4>o 
with once adjacent interchange to be 0. Thus, it is rea¬ 


sonable to define 1{P) as the number of least adjacent in¬ 
terchange from order 4>o to P4>o. This is consistent with 
the field operator language where we are barely allowed 
to interchange neighbor operators one by one. /(P) is 
called length of P in word metric^^ and it equals to the 
number of inversion pairs (nm) in ordered P4>o(anyon-m 
appears on the left of anyon-n in P4>o for m > n)[l^. 
Typically, 1{PP') < 1{P) -f 1{P') and surely is 

not a representation of the permutation group, unless 
e~'^^ = ±1. The negative sign in the phase is to ensure 
the positive interchange phase of the total wavefunction. 
To gain more reliability, we label the anyons with upper 
indices and let P permute on anyons, equivalent to on the 
states, to check the interchange properties of the wave- 
function. /(P) now is the sum over the inversion pairs 
parameter , since those anyons must be interchanged 
to obtain P4>o. For instance, two-anyons 

^ E = 44 + (2) 

P 

and three-anyons 

^llc =444 + e-^0s^-4a44 + e-^<^a--e-^0s^-4a44 

For two-anyons wavefunction, interchange anyon 1 and 2 
denoting as 

PH^4=4a4+e-^‘^^^-44a 

=e-^o^^^i4a4+e^''^^^44) ( 4 ) 

If 921 = —^ 12 , then the interchange yields a phase shift 

~ Certainly interchange twice 

^ ^12 ( 5 ) 

gives no effect. We shall emphasize here that also 

interchange the indices in the phase term. Let gnm = 
—9mn = 1 for n < m in general, corresponding to the 


2 


clockwise/count-clockwise winding. Then the permuta¬ 
tion of anyons leads to a total phase shift (see the proof 
in Appendix. 


For different permutation, the phase shift for the wave- 
function is not identical, otherwise it must be a trivial 
representation of permutation group as Boson. The 
interchange phase of the wavefunction(j 6 ]) is anyons label 
dependent and this is due to the selection of anyons 
order in $o- Different selection of gives different 
wavefunction T' and yields different interchange phases. 
Yet, this does not hurt the indistinguishability and 
statistics since 

In the totally ordered coordinate space Y, the anyon 
state is projected as wavefunction Now 

the coordinate Xi works as label and gij is replaced 
with g{xi,Xj). g{xi,Xj) = -g{xj,Xi) = 1 for Xi ^ Xj. 
Here ^ is the order relation symbol. In the region 
xi ^ X 2 ^ ••• ^ xat, all the involved inversion pair 
parameter g{xi^Xj) = 1 and the wavefunction is 

equation([T]). In other regions we could permute the 
coordinates to increasing order and the wavefunction 
equals to the permutation on wavefunction([T]). Compare 
with the charge-flux model [H, Q, our many-anyons 
wavefunction can ensure the permutation properties 
for many anyons and work for any Hamiltonian system 
manifesting the anyonic interchange phase. 


Capaeity of state-The many-anyons wavefunction ([T]) 
is neither symmetric as Boson nor anti-symmetric as 
Fermion. Consider two and three anyons at the same 
state(a = 6 = c), equations (| 2 ]) and (|3|)becomes 

< i>aa =(1 + e -^ 0 ) cPa < t>a ( 8 ) 

^'ooo =(1 + e *^)(l + e + e ‘^^^)<t>a4’a4’a (9) 

By induction(or see the proof through generating func¬ 
tion in Appendix. H) the wavefunction of Na anyons at 
the same state-a. 


Na — l n 

= n ^ 

n=0 /c=0 

2 _ ^—i*nO 

~ II( 1 - ^ ^ 4’a---4’a 

n=l 


( 10 ) 


Indeed, this is the Q-faetorial form for the number op¬ 
erator in quantum group [Toj with Q = e Anal¬ 
ogy to the exclusion of Fermion, where more than two 
Fermions at the same state vanish the total wavefunc- 
tion(| 8 ]), i.e.l -h = 0 , the many-anyons wavefunc¬ 
tion ([ 10 ]) vanishes when Na ^ ^+ 1 , providing = 

1 . As a consequence, the maximal occupation or eapaeity 


of state-a is q. This is the general exelusion prineiple we 
obtained from the interchange phase! 

In order that a thermodynamic limit can be achieved via 
a sequence of systems with different particle numbers [ 3 , 
[l^ . the phase O/tt shall be a rational number r/p. The 
capacity q of quantum state is always an integer. 


Q = 


2p-l, 
p- 1, 


odd 


r IS even 


( 11 ) 


Irrational phase factor O/ir can not vanish the total wave- 
function, thus yielding no constraint on the occupation. 
In the (real)Q-analogy quantum group approach [Tol- 
21 |(aa’*' + Qa^a = with — 1 < Q < 1 ), it is noticed 
that the number operator is positive definite unless Q 
is the root of unity. If we could extend their Q to the 
anyonic interchange phase e~^^ ^ the capacity of quantum 
state should be found. 

Many-anyons wavefunction with mutual exchange phases 
0o,g among species a, /3 could be similarly achieved. Basi¬ 
cally we construct the wavefunction via equation ([T]) with 
01{P) replaced by 0(P) = {j'l) is inversion pair 

in P^o likewise and ai = a, /3 • • • is the specie anyon-i 
belonging to. For Na anyons at state-a of specie a and 
Nh anyons at state -6 of specie /3, 






N a — In Nij — ln 

n(E n(E 

n=0 k=0 n=0 k=0 


X 


E- 




( 12 ) 


In the last sum of permutation terms, there are 
{Na + N}j)\/(NalNijl) distinct permutations and all the 
mutual exchange phases appear in these terms. 
They are mutual orthogonal and the superposition will 
never cancel each other. Thus the wavefunction can only 
vanish due to the first and second Q-factorial terms of the 
same state interchange phase. Consequently, the occu¬ 
pation of the same species different states((a = /3, a ^h) 
and of the different species (a 7 ^ /3) do not mutually 
affect. Only the interchange phase of the same species 
Oaa constrain the capacity of each single state, quite 
different from the exclusion statistics parameter g^g 
from fractional exclusion statistics [l^. [l^ . 

As to the generalized ideal gas of fractional exclusion 
statistics without mutual statistics [l^ gag = Oidag^ 
specifically, a Fermi-like step distribution is found at 
T = 0. Below the Fermi surface, the occupation number 
is 1/a. The statistics parameter is mapped to the 
phase by a = O/t: = r/p, so the maximal occupation 
ne<EF necessary an integer. This is not the 

same to the capacity we find here either. For instance, 
when 0 = 27r/3, the state capacity is g = 2 while from 
fractional exclusion statistics, = 3/2 and for 

^ = 7 r/ 2 , we get q = 3 while ne<EF — 2- 


Statisties- From our construction of many-anyons 
wavefunction, the explicit statistical effect for anyon is 
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that the interchange phase 0 gives a finite capacity of 
each quantum state. Finite c apacity q of quantum state 
is defined as Gentile statistics |22l - l24| . The grand canon¬ 
ical partition function is 


z=Es:-’.”-=ns: 

{».} * 


i ni=0 


qGj 


= Y, = n E (14) 


{Nj} 


j Nj=0 


(13) 


here Zi = {rii} denote all the possible config¬ 

uration of rii anyons at quantum state-i and {A/’j} de¬ 
note all the possible configuration of Nj anyons at en¬ 
ergy level Cj with degeneracy Gj. The statistical weight 
— T Tj W n and Wj is tl[ms{combination with lim¬ 
ited repetition\2^) 

k=lNj/iq+l)} 

wj= Y ( 15 ) 

k=0 

here Binomial constant Cq = G\/{k\{G — k)\). Then the 
mean occupation is 


{ni)=Y'^iYl 

{rii} r 


Zj (1 + q)z^+‘> 
1-Zi zl+‘‘ - 1 


= {Ni)/Gi 


(16) 


The thermodynamics behavior such as the heat capac¬ 
ity, equation of state and condensation temperature 
are well studied [H, [l^. From equation ([14]), the most 
probable occupation Nj is determinant from the identity 


Sn, log(Ty,-^f) = 6n, logiy^ljv, + log^i = 0 (17) 

Compare equations m and (HZD , we can estimate 
the variance between the mean occupation and the 
most probable occupation A = {Nj)/Nj — 1. It can 
be numerically verified that the variance is almost 
zero (filled square in FIG(T]) coinciding with intuition. 
Certainly, at T = 0, (rii) = q for those states <M, 
similar to the fractional exclusion statistics 13|, |l^ . 
Plausibly, one may consider that Nj particles oc¬ 
cupy at least [Aj/g] (ceiling function) states among 
Gj states, then the effective Fock space dimension 
dB = Gj + Nj - \Nj/q'\ « Gj + {Nj - 1)(1 - 1/g) and 
the number of ways Wj = is the same to Wu’s 

count ingQ, providing q = l/a. 

Yet numerically Wj ^ Wj in general. From identitv ([T7|) . 
the most probable occupation is controled by the 
differentiation on particle number SN^ogW{N). It is 
close for the two different kinds of counting only in 
the limit of Gj ^ Nj as shown in the inset of FIGUl 
This corresponds to the most probable occupation 
number of high energy state and both of them reduce 
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FIG. 1. Numerical comparison of the variance: Ai = 
{N)/N —1 is the variance between mean occupation and most 
probable occupation for Gentile statistics. A 2 = (N)/N' — 1 
compares the mean occupation of Gentile statistics and the 
most probable occupation of fractional exclusion statistics. 
Here we set G = 10^ and g = 3 as an example. The in¬ 
set is the relative difference of particle number differentiation 
S = log lF/(5w log VF' — 1, which represents the difference 
of most probable occupation between Gentile and fractional 
exclusion statistics. 


to the classical Boltzmann distribution. For the low 
energy state, Nj ^ Gj, Gentile statistics and fractional 
exclusion statistics yield completely different occupation 
number (see also empty circle data in FIGHj). 


Lgnon excitation m Heisenberg spin cnam via Hetne 
[iqI . [ 27 1, the charge-flux model 0 and the spin- 
from Jordan-Wigner transformation 


Discussion- It is confusing to talk about many anyons 
at the same state due to the interchange phase, since we 
would have For ^ 7 ^ 0 , it seems that the 

wavefunction vanishes. In some proposed model, such as 
the magnon excitation in Heisenberg spin chain via Bethe 
ansatz 

S on mapping 

people always preinstall the hard-core condition, 
resulting in the conventional Fermi statistics [l^. In¬ 
deed in the coordinate space of N indis¬ 

tinguishable particles, the permutation of coordinates 
at point(xi,X 2 ..Xat) gives the identical point while the 
wavefunction feels a phase shift Q. The wavefunction 
is multivalued [H, Q and the interchange property(T = 
e*^T) does not directly indicate a vanishing wavefunc¬ 
tion. The hard-core condition is a too strong constraint 
for identical particles and it makes sense to discuss many- 
anyons at the same state. 

We also noticed that the capacity of quantum state is 
super sensitive to the interchange phase. Two very close 
phases 0 and 0' might yield completely different capaci¬ 
ties. Thus to obtain a stable capacity and statistics, the 
interchange phase shall be exact. If the phase factor r/p 
could be looked as the ratio of electron and flux, thus 
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the resistance, then the resistance must be exactly on 
the plateau, well known result in FQHEjiSIsil. 

In our construction of many-anyons wavefunction([T]) 
from single anyon state, the permutation is decomposed 
into least adjacent interchange. We applied the winding 
interpretation gnm = —gmn = 1 for n < m and m 
reduces to the length of P in the word metric. Then 
anyons shall obey Gentile statistics. Yet, the permu¬ 
tation property (|28|) is valid for any gnm = —gmn^ e.g. 
gnm = —gmn = (—for Ti < 171. The statistics 
effect may be discussed in future work. 

More generally, the many-anyons wavefunction is 

^'= (18) 

P 

here 4>o could be any function(not necessary to be the 
product of single anyons state) and P permutes on any 
indistinguishable indices. The permutation P' on the 
wavefunction 

p 

p 

^gie'(p') y^g-ie(p'p)p/p^^ ( 19 ) 

P 

if 

0(P'P) - P'0(P) = 0'(P') (20) 

is independent of P, then the permutation P' yield a 
total phase shift 0'(P') for T. Equations (fT8|) and (f2Q|) 
are the general definition for anyons we proposed. There 
might be other form for 0(P) satisfying eg nation (|2Q|) 
and thus controlling other anyonic statistics. 

The many-anyons wavefunction equations ([T]) and ([T8|) 
may be used in the consensed matter system and the 
thermodynamics of Gentile statistics might be compared 
with experiments. 


I.Interchange properties 

For convenience, we write the ordered anyons sequence 
P*-^^^ -4>o. Sequence ijkl.. could represent 
the permutation P. Any permutation could be decom¬ 
posed into least adjacent interchange called reduced word 
in word metric [l^. Although the decomposition is not 
unique, the interchanged anyons pairs are identical. They 
are the inversion pairs (mn) in sequence ijkl..^ if m is 
on the left of n for m > n[l8|. In terms of inversion pair 
parameter gnm 

= Y, 9nm (21) 

{mn) 


Now we will prove the interchange property of the many- 
anyons wavefunction - ' Consider the following two 

sequences(i > j): 

(E)f ~ (l\(l2‘--eLk—l ^ (^k-\-l^k-\-2“(^m — l j 

^^ ^ ^ 


(^2 — (2i(22• •1 j (^k-\-l(^k-\-2"(^m — l ^ (^m-\-l^m-\-2”(^n 

" -V-" "-V-" "-V-" 

si s2 s3 

(23) 


Z(P@i) = E gjmi T ^ ^ gimi ^ ^ gm2j T ^ ^ gim2 

mi>j mi >2 m 2 <j m 2 >i 

+ ^ gmsj + E gmsi + ij-indept terms. 

m3<j m3<i 

(24) 

;(p©2) ^ E gimi T ^ ^ gjmi T ^ ^ gm2i T ^ ^ fo’m2 

mi>i mi>j m 2 <i m 2 >j 

T ^ ^ gm 3 i T ^ ^ gm 3 j T gji T ji-indept. 

m3<i m3<j 

(25) 

ma{(7 = 1,2,3) are the anyons in segment s^ of the se¬ 
quences. Interchange anyon-i and j by T*-^ , sequences 
©2 and @1 are mutually mapped. The phase shift shall 
be equal, 

kj = f^n{p®‘^) - i{p®^) = f^n{p®^) - i{p®‘^) (26) 

The above identity shall be independent of the position 
of i^j and the order of other anyons. The only solution 
for the above equation is gmn = —gnm and the overall 
phase shift is ^ 

m<i 

hj — ^ ^ {gmj T gim) P gij 

m>j 

pP is a 2-circle permutation as an element of the per¬ 
mutation group and it is easy to check 1(TP) = —hj. 
Thus Any permutation could be de¬ 

composed as some independent cyclic permutation P = 
TPk---Ymn---^ Use the same method, we can prove that 
for any P 

p^f^^iei{P)^j ^28) 

Let gnm = —gmn = 1 for n < m as the clockwise/cout- 
clockwise interpretation, then 1{P) is the length of P and 
the wavefunction reduces to equation([T]). Replace 01{P) 
with 0(P) = the many-anyons wavefunction 

with mutual interchange phase can be obtained in 
the same spirit. 
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II.Occupying the same states 


[3] 


For N anyons sequences, denote lN{k) the number 
of permutation sequences with k < C% inverse pairs. 

= N\/{{N—m)\m\) is the binomial constant. Clearly, 
^Ar(O) = 1 and the generating function ^n{x) of lN{k) 
satisfying the recurrence relation [isj 


[4] 

[5] 

[ 6 ] 

[7] 

[ 8 ] 
[9] 


<i!N{,x) = Y,lN{k)x'^ 

(29) 

[10] 

/c=0 

= (1 + a; + H-h n-i{x) 

(30) 

[11] 

[12] 

For A^a-anyons at the same state, equation([T]) 


[13] 



[14] 


(31) 

[15] 

P 


[16] 

= X^iv(0(e“'^)Va...</>a 

1=0 

(32) 

[17] 

is indeed the generating function of lN{k) with argument 

[18] 

[19] 

e = 1, use eQuation(|3Q|). 


Na — l n 

H ^ 4>a.Pa 


[20] 


[21] 

[22] 

n=0 k=0 



“H^ 1 _ g-ie ) ^ 


[23] 

(33) 

[24] 

n=l 


[25] 

For the wavefunction with different states or 

mutual 

[26] 

interchange phases, the identity equation ([12]) 
checked out in the same way. 

can be 

[27] 



[28] 
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